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We present an exact analytical solution of the spectral problem of quasi 

one-dimensional scaling quantum graphs. Strongly stochastic in the classical 
limit, these systems are frequently employed as models of quantum chaos. We 
show that despite their classical stochasticity all scaling quantum graphs are 
explicitly solvable in the form En = f{n), where n is the sequence number of 
the energy level of the quantum graph and / is a known function, which de- 
pends only on the physical and geometrical properties of the quantum graph. 
Our method of solution motivates a new classification scheme for quantum 
graphs: we show that each quantum graph can be uniquely assigned an integer 
m reflecting its level of complexity. We show that a network of taut strings 
with piecewise constant mass density provides an experimentally realizable 

analogue system of scaling quantum graphs. 
05.45.+b,03.65.Sq 

I. INTRODUCTION 

Quantum graphs [1-3] are the "harmonic oscillators" of quantum chaos. Due to their 
structural simplicity they provide a test bed for a large number of properties and hypotheses 
of quantum chaotic systems. Many theoretical investigations, which are difficult to conduct 
for more familiar quantum chaotic systems [4-6], can be carried out explicitly for quantum 
graphs, both in the classical and in the quantum regimes. An example are recently obtained 
spectral formulas [7-10], which provide explicit analytical expressions for the individual 
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quantum energy eigenvalues of a subset of scaling quantum graphs. 

Recently we were able to generalize our methods to the set of all scaling quantum graphs 
[11]. The purpose of this paper is to provide a more detailed discussion and to present 
new results on the spectral statistics and the convergence of our explicit solution formulas. 
We also present a new classification scheme of scahng quantum graphs. We show that it is 
possible to label each scaling quantum graph with an integer m which reflects the degree of 
complexity of its spectrum. We also suggest an experimentally realizable analogue system of 
scaling quantum graphs. This shows that scaling quantum graphs are more than academic 
constructs, and that physical systems can be found which can be analyzed on the basis of 
the theory of scaling quantum graphs. This view is corroborated by a recently published 
microwave realization of quantum graphs [12]. 

Our paper is organized in the following way. In Sec. II we introduce scaling quantum 
graphs and review briefly explicit spectral formulas obtained for a sub-class of scaling quan- 
tum graphs. In Sec. Ill we examine the spectral equation of scaling quantum graphs. In 
Sec. IV we define spectral separators whose knowledge enables the construction of explicit 
spectral formulas for scaling quantum graphs. We also define a new spectral hierarchy of 
scaling quantum graphs which is based on the complexity of their spectra. In Sec. V we 
investigate the spectral statistics of quantum graphs. We show that because of the existence 
of a spectral cut-off the spectral statistics of finite quantum graphs are never exactly Wigne- 
rian. We investigate the spectral statistics of a four- vertex scaling quantum graph in detail. 
Comparing its spectral statistics with the spectral statistics of more highly connected quan- 
tum graphs we show that the index m, although indicative of the complexity of the spectrum 
of a quantum graph, does not uniquely characterize its spectral statistics. In Sec. VI we 
present Lagrange's inversion formula as a new and alternative method for obtaining explicit 
spectral formulas. In Sec. VII we discuss our results. In Sec. VIII we summarize our results 
and conclude the paper. The paper has two appendices. In Appendix A we provide a simple 
proof for the statement that the spectral equation of the m = complexity sub-class of 
scaling quantum graphs has one and only one root per root cell. This is important since 
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our theory of explicit spectral formulas of scaling quantum graphs crucially hinges on this 
statement. In Appendix B we show that our spectral formulas are indeed convergent, and 
in addition that they converge to the correct spectral points. 

II. SCALING QUANTUM GRAPHS 

As illustrated in Fig. 1, quantum graphs consist of a quantum particle moving on a 
one-dimensional network of bonds and vertices. 




FIG. 1. Quantum graph: a quantum particle moves along the bonds of a generic graph and 
scatters at its vertices. 

The bonds Bij of the graph may be equipped with potentials Uij. We refer to these po- 
tentials as bond potentials or the dressing of the graph bonds. The parameters determining 
the strength and the shape of the bond potentials are referred to as dressing parameters. 
In what follows the bond potentials are considered to be scaling potentials, Uij = XijE, 
Xij = Xji = const. The physical meaning and the reason for introducing the scaling assump- 
tion are discussed in [7-10]. In addition, in Sec. VII, we present a physical analogue system 
of scaling quantum graphs, a network of taut strings, which has the same spectral equation 
as scaling quantum graphs. The string system is an example of a naturally scaling system. 
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In a more general context one can consider the scaling assumption as a tool which allows to 
avoid unnecessary mathematical complications. For most physical systems scaling can be 
achieved, even experimentally [13], by an appropriate choice of parameters. We also define 
E — since for the discussion below it is frequently more convenient to work with k than 
to work with E. 

For ^ = quantum graphs produce strongly stochastic (mixing) classical counterparts 
- a classical particle moving on the same one-dimensional network, scattering randomly on 
its vertices [1,2,14-16]. We use the word stochastic to characterize the classical dynamics of 
the particle on the graph since classically the scattering at the vertices is not a deterministic 
process as required for deterministic chaos [17], but a random, stochastic process, where the 
classical scattering probabilities are determined directly from the quantum dynamics in the 
limit h^O [18]. 

Despite the apparent simplicity of quantum graphs, their behavior exhibits many familiar 
features of classically chaotic systems. Examples are the exponential prohferation of classical 
periodic orbits and the approximate Wignerian statistics of nearest-neighbor spacings [1,2] 
(see also Sec. V). As a result quantum graphs are quantum stochastic systems, which mimic 
closely the behavior of quantum chaotic systems. It is therefore very interesting that despite 
their classical stochasticity and despite many familiar phenomenological features of quantum 
chaos exhibited in the quantum regime, the spectral problem for scaling quantum graphs 
turns out to be explicitly solvable [11,19]. 

Let us first outline the solution for a particular class of scaling quantum graphs, called 
regular in [7-10]. We note that the term "regular" as used here refers to the regular behavior 
of the spectrum of the corresponding quantum graphs and has nothing to do with regular 
graphs as defined in graph theory [20], e.g. graphs with a fixed coordination number. A case 
in point is the recent paper by Severini and Tanner [21] where the term "regular quantum 
graphs" refers to quantum graphs with a special graph topology. 

For regular quantum graphs there exists a set of /c-intervals /„, each of which contains 
precisely one momentum eigenvalue kn (see Appendix A). The end points A;„ of these inter- 
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vals, /„ = kn-i, kn , form a periodic set, 

kn^Kin + K2, (1) 

where the constants Ki, K2 are determined exphcitly in terms of the parameters of the quan- 
tum graph. Clearly, the points kn separate the eigenvalues kn from each other, and are 
therefore called separators (see Sec. IV). 

As soon as the separators kn and the density of states p{k) are known, an explicit 
expression for the energy eigenvalues of a given quantum graph is obtained either by first 
computing the momentum eigenvalues 

kn^ p{k)kdk, (2) 

and then using = k^^^ or by computing directly as 

En = r P{E) EdE, (3) 

where En = k^, p{E)dE = p{k)dk. An explicit periodic-orbit expansion of the density of 
states p{k) is given by [7,8] 

p{k) ^Y.^{k-kn)^^ + Re^Yl E ^P^^'^"^'' (4) 

n p v=l 

where S^, and Ap are correspondingly the reduced action lengths and the weight factors 
of the prime periodic orbits labeled by p, v is the multiple traversal index, and Sq is the 
total reduced action length of the graph [9]. The constant term in the expansion (4) of 
p shows that ki in Eq. (1) is given by ki — tt/Sq. In order to illustrate the construction 
of explicit spectral formulas we assume, for simplicity, that K2 — 1/2 and all Ap are real. 
Both assumptions hold for a large class of regular quantum graphs. If we now use the 
expansion (4) in Eq. (3) we arrive at the following exact, explicit periodic-orbit expansion 
of the individual energy levels of the corresponding regular quantum graphs: 



where ujp = t:S^/ Sq. Therefore, according to (5), the index n that counts the separators kn of 
the regular quantum graph, is a quantum number in the sense that it exphcitly enumerates 
the physical eigenstates. In this respect, the explicit formulas for the quantum energy 
levels En of these systems are analogous to the well-known Einstein-Brillouin-Keller (EBK) 
quantization formulas for integrable systems [4-6] . This is a very interesting fact from the 
point of view of the semiclassical periodic-orbit quantization theory. In this respect, the 
regular quantum graphs represent curious hybrids of classical stochasticity and quantum 
spectral solvabihty. 

However, the systems for which the expansion (5) is valid, represent a very special class 
of quantum graphs. Just how special such "spectral regularity" is, can be illustrated in 
terms of the behavior of the corresponding spectral staircase function, 

N{E)^Y.^{E-En), (6) 

n 

where 9 is the unit step function defined as 

{0, for X < 0, 
1/2, fora; = 0, (7) 
1, for X > 0. 

It was shown in [7,10], that for the regular systems, the average spectral staircase (Weyl's 

average), 

7V(^) = — /E + iV(0), (8) 

TT 

has the piercing property, i.e. it intersects every stair step of the spectral staircase function 
N{E), as illustrated in Fig. 2. 
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FIG. 2. Piercing property of the regular quantum graphs. The spectral staircase function of a 
regular quantum graph is pierced by its average N{E). 

If a quantum system has the piercing property, there exists exactly one intersection point 
En-i, between every two neighboring energy levels En~i < En, 

N{En-i) = N{K^i), K-i < K-i < n = 1, 2, . . . . (9) 

The En thus defined may serve as separators for the quantum energy spectrum. As shown in 
Fig. 2 the piercing- average requirement (9) is indeed quite restrictive. Consequently, regular 
quantum graphs form a relatively small subset of quantum graphs. As demonstrated in 
[10,22], only a few graph topologies (for instance linear chains) admit a regular regime 
for an appropriate choice of network parameters. As an example, a four-vertex linear-chain 
quantum graph (see inset of Fig. 3), which is characterized by the values of the two reflection 
coefficients r2 and at the two middle vertices V2 and V3, is in the regular regime if these 
parameters fall into the shaded region shown in Fig. 3. 
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The majority of scaling quantum graphs do not admit regular regimes. Hence it is 
intriguing to understand the spectral behavior of irregular quantum graphs, i.e. those for 
which the piercing-average condition (9) is violated. 



III. SPECTRAL EQUATION 

In order to set the stage for the following discussion, let us recall some general definitions 
and properties of quantum graphs. As mentioned in the introduction, a quantum graph [1-3] 
consists of a quantum particle moving on a one-dimensional network of Nb bonds connecting 
Ny vertices (Fig. 1). Every bond Bij which connects the vertices Vi and Vj, carries a solution 
of the Schrodinger equation, Hipn = Eni'n- The length of the bonds is denoted by Lij. With 
the constant scaled potentials Uij = XijE defined on the bonds of the graph, the Schrodinger 
equation is 

^^Pi,{x)+(3lE^P,,{x)=0, (10) 

where Pij = ±a/1 — Xij. 

Below we shall assume for simplicity that the energy E is kept above the maximal scaled 
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potential height, i.e. Xij < 1, i,j — 1, ...,Nv, so that tunnehng solutions are excluded and 
the general solution of Eq. (10) on the bond B^j is 

ilJij{x) = aijC'''"^'^'''^ + bijc''"^*^'''^ . (11) 

The quantization conditions for quantum graphs are the result of the requirement that the 
solutions (11) must satisfy the continuity and the current conservation conditions at every 
vertex V^. The procedure of imposing the boundary conditions can be reformulated in terms 
of an auxiliary problem of quantum scattering on the vertices of the graph [2,10,14], which 
provides an elegant solution of the graph quantization problem. As shown in [2,10,14] the 
consistency of the complete set of boundary conditions at all vertices yields the spectral 
equation 

A{k) = det [1 - S{k)] = 0, (12) 

where S{k) is a 2Nb x 2Nb unitary (scattering) matrix [2,10,14], 

Sij{k) = tije'^'^^'K (13) 

Here the capital indices /, J are used to denote the directed bonds, I,J— 1,2,..., 2Nb. We 
denote by /' the time-reversed bond of /. The elements tu (discussed in detail in [10]) have 
the meaning of transmission (reflection) amplitudes for transitions between the (directed) 
bonds / and J. Transmission occurs if / and J are connected and J ^ I' . II I and J are 
not connected, we have tu = 0. An example here is tu = for all /. For J = I' the 
matrix element tu has the meaning of a reflection amplitude [2,9,10,14]. Due to the scahng 
condition, the i/j's are constant (fc-independent) parameters. 

For conventional quantum graphs without potential dressing the connection between the 
coefficients tjj and the expansion coefficients Ap in Eqs. (4) and (5) was established early 
on in the seminal literature on quantum graphs, e.g. in Refs. [1,2]. Later it was shown to 
hold also in the case of dressed, scaling quantum graphs [15]. Each transition of an orbit p 
from a bond I to J contributes the factor tu to the weight Ap of the orbit, so that 
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Ap^lltij, (14) 
{p} 

where the product is taken over the sequence of bonds traced. 

Note that the phases of the exponentials in Eq. (13) coincide with the classical actions 
associated with the particle path traversing the bond Bj, 

Si(k) = PiLjk. (15) 

The spectral determinant (12) can be written in the form 

A(fc) = e^®°«A«(A;), (16) 

where /S.ji{k) is the (real) modulus of A(A;) and ©o(^) is its phase. The phase is given by [9] 

Qoik) = ^ In det 5 = Sok - n^o, (17) 



where 5*0, the total reduced action length as introduced in Eq. (4), is given exphcitly by 

2A'b 
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So^-Y^LiPj (18) 



7=1 



and 7o is a constant phase. The modulus is given by [9] 

Nr 

Aij(/c) = cos {Sok - 7r7o) - ^ cos{Sik - n-fi), (19) 

1=1 

where Oj are constant coefficients, ji are constant phases, A^r is the number of harmonic 
terms in the sum of Eq. (19) and the frequencies Si are hnear combinations of the reduced 
classical bond action lengths Sj — PiLj. Sq is the largest frequency in Eq. (19), i.e. Si < Sq, 
i — 1 ... Nr. This fact will be of crucial importance below. 

The spectrum of the quantum graph is obtained from the equation 

Anik) = 0. (20) 

In Appendix A we prove that if the coefficients of the characteristic function $(/c) of the 
graph, 
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$(A;) = y^g^ cos{Sik 



(21) 



i=l 



satisfy the condition 



\ai\ = a < 1, 

i=l 



(22) 



precisely one solution /c„ of Eq. (20) can be found between each two sequential separators 



fc„ = (n + 7o + + 1) 

^0 



(23) 



where /x, an integer, is to be adjusted such that ki < ki < k2. This is the case, e.g. for a 
two-bond graph (Fig. 4) with the bond lengths Li and L2, for which the spectral equation 



IS 



sin(5'o/c) — rsm{Sik) = 0. 



(24) 



Here Sq = LiPi + L2P2, Si = LiPi — L2P2, and r is a constant positive reflection coefficient 
at the vertex V2 between the two bonds. Since |r| < 1, the condition (22) is satisfied and 
hence this graph is always regular. 

20 




Momentum, k 

FIG. 4. The three-vertex linear graph (inset) and the corresponding staircase function. The 
intersections between N{k) and N{k) correspond to the separating points kn- 
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In this case every step of the spectral staircase function (6) is pierced by its average 
(Fig. 4), or equivalently, every interval J„ = [kn-i, kn] contains precisely one quantum eigen- 
value of the momentum. This spectral regularity is the key for obtaining the explicit har- 
monic expansion for each individual root of the spectral determinant (12). In general, 
however, the regularity condition (22) does not hold and hence the principle "one root per 
interval J„" (see Appendix A) is violated. This is illustrated in Fig. 5, which shows the 
behavior of the spectral staircase for the four-vertex linear chain in two different dynamical 
regimes. The spectral staircase on the right corresponds to a case in which the parameters 
r2 and fall outside of the shaded regularity region in Fig. 3. 
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FIG. 5. The staircase and the average N{k) for the four-vertex hnear graph in a regular (left) 
and in an irregular (right) regime. 

Hence, in order to proceed with an analysis similar to the one for regular quantum graphs, 
one needs to find a set of separating points that "bootstrap" the spectrum, and allow us to 
integrate around each delta-peak of p{k), as in Eq. (2). 

IV. SEPARATORS 

What is the set of points that can be used as separators for a generic quantum graph? 
Since the points kn that need to be separated, are the zeros of the spectral determinant 
Aji{k), one can invoke an elementary, classic theorem [23], which states that between every 
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two roots of a real, continuous, different iable function there exists an extremum point. 
Moreover, extending A/j(fc) into the complex plane, Ar(/c) ^{z), and using the Hadamard 
representation of the resulting entire function ^{z), 



A(z) 



n 1 



(25) 



where g > is the multiplicity of the root = 0, and all the roots kn are assumed to be 
real as required, since A/j(A;) is derived from a Hermitian eigenvalue problem, one can show 
[24,25] that there is exactly one zero between every two neighboring extrema of A/j(/c), i.e. 
that the zeros and the extrema interlace and "extra wiggles" such as, e.g., illustrated by the 
dashed line in Fig. 6, are not possible. Hence the locations of extrema may be used as the 
separating points for bootstrapping the physical spectrum. 




FIG. 6. The interlacing sequence of roots and extrema of the spectral determinant. The dashed 
line represents the forbidden "extra wiggles". 

Strictly speaking, all this works only for simple roots of A(2;), which is the generic case. 
Multiple roots may, and in special cases do, occur. But as explained in Refs. [11,19], these 
cases are trivial to deal with. In such a case the separators and the spectral points are 
degenerate and no further computation is necessary. 

Is it any easier to obtain the extrema of /S.fi{k) than to obtain its roots? Interestingly, 
looking for the answer to this question provides us with a complete scheme for establishing 
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a hierarchy of quantum graphs according to their spectral irregularity. Let us examine this 
question more closely. 

The equation A^(A;) = that defines the extrema of Ar(A;) is 

Nr 

sin (^oA; - ttjq) - ^ a^ei sm{Sik - irji) = 0, (26) 

i=l 

where 

= f < 1- (27) 

Note that this is the same type of trigonometric polynomial as the original spectral equation 
(20) (with a shifted phase 70), with the new characteristic function 

Nr 

$(1) (k) = J2 <^i^i MSik - tttO . (28) 

i=l 

However, compared to the original $(/c) of Eq. (21), this function has certainly a better 
chance of satisfying the regularity condition (22), because the absolute values of the coeffi- 
cients Gi have been decreased by the factors < 1, i.e. 

Qi aiCi. (29) 

Let us assume that for a certain irregular graph with ^1^=1 > the new characteristic 
function ^^^\k) actually does satisfy the regularity condition (22), i.e. 

J2M<'^- (30) 

i=l 

According to the results of Sees. II and III this means that the zeros of A^(/i;) (the extremal 
points of Aji{k)) can be separated from each other by the periodic sequence of separators 
(23), i.e. that there exists exactly one extremum of Aij(/c) between every two points kn- For 
functions of the type (25) the converse statement is also true [24,25], i.e. there exists a root 
of A.ji(k) between every two extrema of Ar(/c). This suggests a direct strategy for obtaining 
the roots of Eq. (20). First, as mentioned above, we note that the spectral equation for 
A'j^{k) can be written in the form 
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Ny 

Kik) = E + '^^•) = 0' (31) 

i=o 

where Cj and (pj are constants. From this we obtain the following explicit formula for the 
density- of- extremas functional p^^\k) 



pW(^) = iAMi5(Ai,(^)) = iAMi- y 



oo 



oo oo „oo -^r 



no=— oo riiVp 



JV- 



r 



exp^i^ni[5iA; + (^i] L (32) 



1=0 



where the Jn are Bessel functions of the first kind [26], and the integrals in Eq. (32) converge 
for Cj 7^ 0. Now, using the expansion (32) together with the periodic separators (23), one 
obtains the separating points kn^ for the roots A;„ of Aii{k), via 

L P^'\k)kdk. (33) 

Here we used the notation fcn'* for the periodic separators (23), for future convenience. 
Following this step, using the separators kn^ obtained in Eq. (33), we find the roots A;„ of 
the spectral equation via 



£(0) 



k^= / p^''>(k)kdk, (34) 

where the notation p^^\k) was used for the density of states p{k). 

For the case of the four-vertex linear graph, this situation is again illustrated in Fig. 3, 
in which it is now assumed that the differentiated equation (26) satisfies the regularity 
condition all through the domain — 1 < r2, rs < 1. This would be the case, e.g., if the bond 
action lengths are chosen to be 5"^ = 0.25, 5"^ 0.45, and = 1 - S'} - S^. For this case 
one can immediately verify that the spectral equation of the four-vertex linear chain, 

sm{Sok) = rs sm{S^k + S^k - S^k) - r^r^ sin(5°A; - S^k + S^k) 

+r2 sivL{Slk - Slk -S^k), (35) 

is irregular outside of the region |r3| + |r2r3| + |r2| < 1, but the coefficients of the differentiated 
equation satisfy the regularity condition (22). 
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Clearly this strategy can be applied in the general case. If the differentiated equation 
(26) is not of the regular type, one can differentiate the spectral equation (20) as many times 
as it is necessary to obtain an equation of regular type at the m-th step. Indeed, the l-th 
derivative of the spectral equation is 

JVr 



A^S{k) = cos ^^0^ - TTTo + - ^ Oje- COS (^Sik - W-fi + 



(36) 



Obviously, since all e^'s are smaller than 1, we eventually (after a finite number m of steps) 
arrive at an equation that satisfies the regularity condition (22), 

5^|a,en<l. (37) 

1=1 

An upper bound for m is easily established. We have m < — ln(^^-^ |aj|)/ ln(maxj e^). 
Then, once the condition (37) for the m-th derivative of its spectral determinant is satisfied, 
its zeros are separated from each other by a periodic sequence of points, 

^i™^ = ^(n + 7o + /x+l) (38) 

•JO 

as m Eq. (23). Using the density p^'^^k) of zeros of A^^\k), which is obtained explicitly in 
complete analogy with Eq. (32), we can evaluate the zeros themselves as 

U^-^^ = / p^"'\k)kdk. (39) 

Obviously, these points are now the extrema of and moreover, since there is 

exactly one root of A^^(A;) between any two neighboring points A^""^^(/c) has no 

extrema other than the ones obtained in Eq. (39). 

The newly obtained extrema of A^~^-*(A;) separate its zeros from each other, and hence 
serve as the separators for the roots of A^~^^(A;). As a consequence, we can now find all the 
roots of A^~^'*(A;) by using Eq. (39) recursively until we arrive at the 0-th level to obtain 
the roots A;„ of the original spectral determinant. This solves the problem of obtaining the 
energy spectrum of any scaling quantum graph. 
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It is important to realize that a quantum graph of a certain topology can still have 
different degrees of irregularity depending on the network's bond lengths and dressing pa- 
rameters (for instance the values of the bond potentials). This point is easily illustrated by 
once more using the example of the four- vertex linear chain. Although the regularity region 
for this graph is always the same, the surrounding blank region in Fig. 3, which corresponds 
to the irregular regime, now acquires structure. Figure 7 shows the spectral diagram for the 
four-vertex linear-chain graph for two different sets of graph parameters, corresponding to 
two different irregularity regimes. 




'2 

FIG. 7. The spectral regime diagram of the four-vertex linear-chain quantum graph. The 
bond-action lengths are chosen to be S\ = 0.2, S2 = 0.6565, 5*3 = 1 — Si — S2 (left panel), 
and Si = 0.1, 5*2 = 0.8565, ^3 = 1 — S'j' — ^2 (right panel) resulting in a maximum degree of 
irregularity of m = 2 (left panel) and m = 6 (right panel). 

The central diamond-shaped regions (m = 0) in Fig. 7 are the same as in Fig. 3 and 
correspond to the same regularity region as in Fig. 3. The outer layers of the regular region 
in Fig. 7 correspond to parameter values that guarantee first (m = 1), second (m = 2), . . . , 
degree of irregularity. 
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V. SPECTRAL STATISTICS 



It is well known that the statistical properties of the spectra of generic quantum graphs 
are well described by Random Matrix Theory (RMT) [1,2]. In particular, the numerically 
obtained nearest-neighbor distribution P{s) of the normalized spacings s [5] of the eigenval- 
ues of highly connected quantum graphs follows closely the profile of the Gaussian random 
matrix ensembles, both in the presence of the time reversal symmetry (GOE), where the 
nearest-neighbor distribution is given by [5] 



and in the absence of it (GUE), where the nearest- neighbor distribution is given by 



This circumstance is one of the most important motivations for studying quantum graphs 
in the context of quantum chaos theory. In particular it is hoped to gain more insight 
into the connection between chaos and random matrix theory and, if possible, to prove the 
Bohigas-Giannoni-Schmit conjecture [27,28], which states that, generically, the spectrum of 
quantum Hamiltonian systems chaotic in the classical limit should conform with the spectral 
properties of the random matrix ensembles. 

However, the exact results presented above show that for quantum graphs with a finite 
number of bonds and vertices (finite quantum graphs) the correspondence with the nearest- 
neighbor distributions (40) and (41), respectively, can only be approximate. Indeed, the 
existence of root separators implies that the eigenvalues of the momentum, kn, will always be 
confined within the root cells, kn G [kn-i, kn\- Hence, for finite quantum graphs, even though 
they may be highly connected, the statistical distribution P(s) of the nearest neighbor 
separations, s„ = A;„ — will be restricted to the finite domain < s < s^ax, and 

will not have the characteristic long tail of the nearest- neighbor distributions (40) and (41), 
respectively. This general property of the spectra of scaling quantum graphs also follows 




(40) 



^gue(s) = exp(-4sV7r). 



(41) 
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from the fact that their spectral function, Afi{k), is an almost periodic function of the 
momentum, and hence its zeros form an almost periodic set [29]. It is clear, therefore, that 
the distances between neighboring points of this set are bounded, i.e. s < Smax, and P{s) is 
zero for s > Smax- These observations, of course, do not preclude the possibility that certain 
finite, highly connected quantum graphs are well, or indeed even exactly described by the 
finite matrix ensembles [5]. 

However, the higher the degree of irregularity of a quantum graph, the larger Smax- Going 
upwards in the "hierarchy of separators" leads to an increase in the allowed nearest neighbor 
spacings, since the maximal possible distance between neighboring separators grows by one 
unit of mean spacing when going from complexity level m — 1 to complexity level m. The 
mechanism for the increase of the allowed maximal nearest-neighbors spacing as a function 
of m is illustrated in Fig. 8. 

*^ • — A'"~''kw^-c\ 



So 



5^ 
^ So 



So 



So 



^ a""\.) = o 



So 

FIG. 8. Increase of the maximally allowed nearest-neighbor separation as a function of m, gen- 
erated by the hierarchy of the separators. At the regular level (A^™"-*), the periodically spaced 
separators kl^^^ (marked by x) separate the roots of A^''(A;i™''') = (full circles). Their maximally 
allowed distance is 2^. The second order separators (the roots of ^\kn'^ = 0), may 
be maximally as far as 3^ apart, etc. The higher the hierarchy of the separator, the larger the 
maximally allowed spacing of nearest neighbors. 



Figure 8 also shows that the roots of a spectral equation A/{(A;) = with irregularity 
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degree m, may be no more than (m + 1)t:/ Sq apart. This provides a simple rule for finding 
an upper limit for Smaxi 

41<rftl = f (^ + 1)- (42) 

Clearly, the possibility of having large separations between the nearest neighbors is necessary 
for producing a statistical distribution for s„ = — kn-i that resembles a Wignerian 
distribution profile, similar to the ones which were numerically obtained in Ref . [2] . 

On the other hand, it is essential to realize that a high irregularity degree m is not enough 
to guarantee Wignerian-like statistics of the nearest neighbor spacings. A simple numerical 
experiment with the spectral equation (35) shows that the separations between nearest 
neighbors do not necessarily assume the largest possible values (42). Hence the degree of 
irregularity indeed provides only an upper limit for the nearest-neighbor separations, and 
does not determine by itself their actual values. 

For example, the dressing parameters of a quantum network can be changed continuously 
so that the system undergoes a transition from an irregularity m regime to an irregularity 
m + 1 regime. As this transition happens, the roots of the spectral equation do not re- 
spond to produce an abrupt increase of the nearest- neighbor separations by tt/ Sq. Instead, 
the maximal nearest-neighbor separation increases smoothly as a function of the dressing 
parameters. 

There is a convenient way to illustrate this increase for the four-vertex chain network, 
using the structure of its spectral regime diagram (Fig. 7). As shown in Fig. 7, the parameter 
regions that correspond to different irregularity degrees for this graph form a system of 
nested diamond shapes, with high irregularity regimes concentrating toward the corners of 
the diagram. A specific set of the action length values, S2, S^, defines the frequencies 
in (35) and hence the maximal irregularity degree rrimax, i-e. the total number of diamond- 
shaped regions, while a choice of the refiection coefficients, r2 and r^, puts the system onto 
a particular point in the diagram. Hence, one can study the effect of increasing irregularity 
by traversing the spectral regime diagram from its center (r2 — — 0) to one of the corners 
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(say, r2 = rs = 1) along the line r2 = = r, < r < 1. For each value of r = r2 = that 
corresponds to a particular irregularity degree, m, one can obtain numerically the maximal 
separation distance, Smax, between the nearest neighbors, and then follow its change as m 
increases. 

In addition to the maximal separation Smax there also exists a minimal separation Smin- 
The vertical bars in Fig. 9 represent the possible range of nearest-neighbor spacings Smin < 
s < Smax for givcu m. Clearly, the maximal root separation is increasing with growing m. 
However the increase is slower than the one given by the linear estimate d^ax in Eq. (42). 
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FIG. 9. Range of nearest-neighbor spacings s/tt as a function of the irregularity degree m for 
the four-vertex chain. The maximal separation was obtained based on the roots found in the 
interval < < vrlO, 000 in each m regime. The bond action lengths are = 0.1, = 0.8999, 
5*3 = 0.0001, which produce a maximal irregularity degree of fTimax — 

27. 

Since the spectral equation (35) is an almost periodic function, the maximal root sep- 
aration found on a sufficiently large finite interval of the momenta (large compared to the 
smallest almost-period of the function (35)) is indeed the maximal root separation produced 
by this function on arbitrary intervals. 

It is also important to notice that the maximal nearest-neighbor separations Smax can 
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be different for two graphs with the same degree of irregularity. Moreover, two quantum 
graphs with the same irregularity degree may have completely different spectral statistics. 
This can be seen from comparing the cases of the topologically simple four-vertex chain 
graph with the fully connected four-vertex quadrangle. The spectral statistics provided by 
the latter example were previously discussed in Ref. [2]. It was shown that the nearest- 
neighbor distribution follows quite closely the anticipated Wignerian shape (both in the 
GOE and in the GUE cases). 

The analysis of the spectral equation of a four- vertex quadrangle with no bond potentials 
and comparable bond lengths (the case considered in Ref. [2]) produces an irregularity degree 
that usually does not exceed m ^ 25. This level of irregularity can be easily achieved by 
the four-vertex chain, which, unlike the quadrangle, does not produce the characteristic 
Wignerian distribution profile. The distribution produced by a four-vertex chain in different 
irregularity regimes is shown in Fig. 10. 
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FIG. 10. Nearest-neighbor spacing probability distributions for the eigenvalues of the four-vertex 
linear-chain quantum graph with bond action lengths Si = 0.1, ^2 = 0.8999, ^3 = 0.0001, in 
different regimes of irregularity m as a function of spacing s (in units of tt). The distribution 
profile is not Wignerian. The higher the irregularity index m, the higher the peak of the distribution 
around s/ir w 1 /5'2 . 
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Some general features of these distribution curves can be easily explained with the help 
of elementary quantum-mechanical arguments applied to the four-vertex chain. Indeed, it 
is clear from Fig. 7, that high irregularity degree for a four-vertex chain can be achieved by 
selecting both reflection coefficients \r2\ and jraj close to 1. Physically, such a choice implies 
that the bonds of the chain are essentially isolated, since the particle almost never transmits 
from one bond to another. This "bond decoupling" also manifests itself in the spectral 
properties of the system by the emergence of three apparent sub-sequences of eigenvalues, 
each associated with one of the "isolated-bond spectra", nn/Si. Hence, one would expect 
that in the case \r2\, Ir^] ~ 1 the nearest-neighbor separations will mostly concentrate around 
the values determined by the inverse bond lengths, s ~ \{7rni/Si) — {irrij/Sj)], where rii, rij 
are independent integers, rather than around a peak defined by the Wignerian distribution. 

Overall, the results of the statistical analysis of the four- vertex chain spectrum show that 
both the small-s and the large-s ends of the P{s) distribution profile change slowly with 
increasing irregularity degree. Even in the case of high irregularity degree, the behavior of 
the roots of Eq. (35) is too restricted by the simple analytical nature of Eq. (35) to exploit 
the possibility of getting as close to, or as far from, one another as is allowed by the hierarchy 
of the separators. 

Since the irregularity hierarchy presented in previous sections is a completely general 
structure, the irregularity degree m produced by this scheme is a very general index. How- 
ever it does not, by itself, determine the spectral characteristics of a given quantum graph. 
While a small irregularity degree can be provided only by a few classes of graphs with rel- 
atively simple geometry, a large degree of irregularity can be shared by a wide variety of 
graphs, which include both the topologically simple ones (with appropriate dressings) and 
the topologically elaborate networks. It is natural, therefore, to expect that the statistical 
spectral properties produced by topologically simple graphs can differ from the ones pro- 
duced by topologically complex networks, even if they are characterized by the same degree 
of irregularity m in the sense of the bootstrapping scheme presented above. 

The four-vertex chain graph, whose spectral equation (35) contains only four oscillat- 
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ing terms, is certainly too simple to produce random-matrix-like behavior, whereas a four- 
quadrangle, whose spectral equation written in the form (35) contains about 830 terms, 
is already sufficiently complex. This situation emphasizes the fact that the general phe- 
nomenological statement "classical chaos implies Wignerian statistics", implicitly assumes 
sufficient complexity of the underlying classical system. 

From the opposite perspective, it may be considered a curiosity that simple networks, 
such as the four-vertex chain, are capable of producing highly irregular spectra. It is inter- 
esting in this context to look for a more refined scheme that could distinguish between the 
complexity of the spectra provided by simple graphs (e.g. linear chains) and the spectra 
of more complicated networks, which are capable of producing random-matrix-like spectral 
statistics. 

In this section we studied spectral properties of quantum graphs only as far as relevant 
in connection with our new "m-scheme". Much more is known about the spectral prop- 
erties of quantum graphs in general (see, e.g., [1-3,21,30-36]). In particular the thrust in 
the investigation of spectral properties nowadays is on understanding spectral correlation 
functions [1-3,33,34,36] and even deriving explicit formuls for them [30,35]. 



The periodic orbit expansions presented in Sec. IV are not the only way to obtain the 
spectrum of regular quantum graphs explicitly. Lagrange's inversion formula [37] offers an 
alternative route. Given an implicit equation of the form 



Lagrange's inversion formula determines a root x* of Eq. (43) according to the explicit series 
expansion 



VI. LAGRANGE'S INVERSION FORMULA 



X — a-\- w(fi{x), 



(43) 




oo 



(44) 



provided (f{x) is analytic in an open interval / containing x* and 
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\w\ < 



X — a 
ip{x) 



y xel. (45) 



Since the regularity condition (22) ensures that the condition (45) is satisfied, we can use 
Lagrange's inversion formula (44) to compute explicit solutions of regular quantum graphs. 

In order to illustrate Lagrange's inversion formula we will now apply it to the solution 
of Eq. (24). Defining x — Sok, the nth root of Eq. (24) satisfies the imphcit equation 

Xn — Tm + (— 1)" arcsin[r sin(pa;„)], (46) 

where p = Si/So and \p\ < 1. Choosing So = 0.3 + 0.7/\/2, Si = 0.3 - 0.7/\^ and 
r = (V2-1)/(V2 + 1), we obtain 0;^'"^''^ = 3.26507..., ar^o"'''*^ - 3L24664 . . . and x^qo"'^ = 

313.98697 We now re-compute these values using the first two terms in the expansion 

(44). For our example they are given by 

42 ) = Trn + arcsin[r sin(p7rn)] ( (-1)" + ^^^"^(pTrn) ^ 1 ^^^^ 

[ yl — sin^(p7m) J 

We obtain a;f ^ = 3.26502 . . ., = 31.24650 ... and xf^^ = 313.98681 . . ., in very good 

, .,T (exact) (exact) i (exact) 

agreement with x\ , x\q and x\qq . 

Although both Eq. (5) and Eq. (44) are exact, and, judging from our example, Eq. (44) 
appears to converge very quickly, the main difference between Eq. (5) and Eq. (44) is that 
no physical insight is gained from Eq. (44), whereas Eq. (5) is tightly connected with the 
classical mechanics of the graph system providing, in the spirit of Feynman's path integrals, 
an intuitively clear picture of the physical processes in terms of a superposition of amplitudes 
associated with classical periodic orbits. 



VII. DISCUSSION 

The first announcement of explicit periodic-orbit expansions of the spectrum of regular 
quantum graphs [8] was universally met with disbelief and puzzlement. It seemed impossible 
to obtain explicit solutions for a quantum system that had been shown to be an excellent 
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model of quantum chaos [38] and, moreover, is completely stochastic in its classical limit [1] . 
However, we found that the rejection of our results was almost always based on the common 
misconception of the "unsolvability" of chaotic systems. We point out here that it is not 
true that classically chaotic systems are necessarily unsolvable. We hope that this insight 
will eliminate much of the reservations commonly expressed toward our results. 

Examples of explicitly solvable chaotic systems are readily available. The shift map 
[39,40], 

= (2x„) modi, x„ e R, n = 0,1,2,..., (48) 

for instance, is "Bernoulli" [40], the strongest form of chaos. Nevertheless the shift map is 
readily solved explicitly, 

Xn^{2''xo) modi, e R, n = 0,1,2,.... (49) 

Another example is provided by the logistic mapping 

Xn+i^ fJ^Xnil- Xn), x„e[0, 1], < ^ < 4, n = 0, 1, 2, . . . , (50) 

widely used in population dynamics [39-41]. For n — 4 this mapping is equivalent with the 
shift map [42] and therefore completely chaotic. Yet an explicit solution, valid at = 4, is 
given by [42]: 

Xji — sin^ (2" arcsin y^) , xq G [0, 1]. (51) 

Therefore, as far as classical chaos is concerned, there is no basis for the belief that classically 
chaotic systems do not allow for explicit analytical solutions. Our contribution in this paper 
is to show that scaling quantum graphs provide the first examples of explicitly solvable 
quantum stochastic systems. 

In this paper we focussed on scaling quantum graphs mainly because of their mathe- 
matical simplicity. However, we will show now that for some physical systems the scaling 
property arises naturally as a consequence of the underlying physics. 
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Consider a taut string of length L, clamped at both ends, with a piecewise constant mass 
density fi{x) = e{x)fio, e{x) = e,, Xi^i < x < Oj, i = 1, . . . , 4, oq = 0, 04 = L, where /Xq is the 
average mass density of the string. This system contains the same physics as a four-vertex 
linear scaling quantum graph since the transverse acoustic excitations of the string satisfy 
the same spectral equation as a four-vertex linear scaling quantum graph. The reason is the 
following. For small transverse oscillations the string obeys the wave equation 

i{x) = 0, (52) 

where ^{x) is the amplitude of the transverse acoustic field of the string at point x and T is 
the tension in the string. 

Equation (52), supplemented with the boundary condition ^(0) = ^(L) = 0, can be 
written in the form (10) of a four-vertex scaling quantum graph. Defining E = lu'^hq/T, we 
obtain 

l^ + A^^l A{xi)^0, (53) 

where /3j = y/el, xi e [ai_i,ai], ipi{xi) = ^(x), a^-i < x < ai. It is obvious that a web of 
taut strings with more complex connectivity as in our example is capable of simulating any 
scaling quantum graph. 

Although the string model has not yet been realized experimentally, a different model has 
been implemented recently in the laboratory [12]. This experiment models a quantum graph 
with the help of interconnected microwave wave guides. The experimental conditions are 
arranged such that only the TEM mode [43] can propagate in a frequency range from about 
100 kHz to 16 GHz. This allows the authors of Ref. [12] to study the spectral properties 
of these microwave graphs in great detail. Even the time-reversal violating case is realized 
with the help of Faraday isolators [44,45]. We suggest here that the authors of Ref. [12] 
could easily modify their experimental set-up to include the case of scaling quantum graphs 
in their measurements. This is done by filling the coaxial cables representing the edges of 
the quantum graphs with dielectrics of different dielectric constants e, respectively. This 
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simple modification would allow the authors of Ref. [12] to extend the set of experimentally 
accessible wave graphs enormously. In addition to the analogues of "conventional" quantum 
graphs (quantum graphs without additional potentials on the graph edges), they would also 
be able ot study the spectral characteristics and periodic-orbit structure of general scaling 
microwave graphs, which are the analogues of scaling quantum graphs. 

The paper by Berkolaiko and Keating [30] is relevant in the context of arriving at explicit 
formulas for physical and mathematical characteristics of quantum graphs. Berkolaiko and 
Keating's result [30], however, pertains to arriving at an explicit formula for the spectral 
form factor K{t) [2,30], whereas the central result of our paper is to present explicit formulas 
for the spectrum itself. In addition the results of Berkolaiko and Keating are derived for the 
special case of conventional, undressed star graphs, whereas our formulas hold for a more 
general class of dressed quantum graphs without restriction of the graph topology. Therefore 
the methods and the physical quantities computed in Ref. [30] are fundamentally different 
from the methods and physical quantities computed in our paper. This also gives us the 
opportunity to clarify a common confusion. It has been suggested to us that our method 
of separators is the same as the method of partitions used in the paper by Berkolaiko and 
Keating [30] , when in fact these two methods have nothing in common. Our separators are 
real numbers which isolate spectral points. The partitions used by Berkolaiko and Keating 
are combinatorial entities related to the number of ways one can represent an integer as 
a sum of other integers. Partitions are a highly interesting mathematical subject, and the 
greatest mathematicians, including the famous Indian mathematician Ramanujan [46] have 
proved deep theorems about them. However, it is clear that both methods are completely 
different, since even from the outset the mathematical categories of the quantities involved 
are different. 

Of particular importance for our investigations is the paper by Barra and Gaspard [35]. 
These authors arrive at an explicit formula for the nearest-neighbor spacing distribution P{s) 
of quantum graphs. Even more. Since the methods of Barra and Gaspard are only based 
on the quasi-periodicity of the spectral equation, their results apply to all quantum systems 
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with a quasi-periodic spectrum, for instance to the dressed quantum graphs discussed in this 
paper. Since our methods yield exphcit formulas for the spectral eigenvalues themselves, we 
hope to be able, in future work, to present alternative explicit representations of P{s) based 
on our explicit periodic-orbit expansions of the spectrum. 

The standard tool of the semiclassical theory used for studying quantum chaotic spectra 
is the periodic orbit expansion for the density of states. Using the density of states approach, 
the individual energy levels are obtained indirectly, typically with semiclassical accuracy, as 
the singularities of the periodic orbit sum. For quantum graphs, however, it turns out that 
one can go one step further, and express the individual quantum energy levels En in terms 
of exact, explicit formulas. Moreover, energy levels can be targeted and labeled individually 
and computed individually without the necessity of knowing any of the preceding energy 
levels. In addition we showed that we can assign a unique degree m to any given quantum 
graph, where m defines the minimum number of differentiations of the spectral determinant 
necessary to reach the regular level, which bootstraps the spectrum. Thus quantum graphs 
appear to have a certain intrinsic degree of complexity which is characterized by m. 

As discussed in Ref. [10], in order to obtain the expansion (5) for a generic quantum 
graph, one needs to obtain the piercing average of the spectral staircase, which, in general, 
is a complicated task. The proposed scheme for bootstrapping the spectrum represents a 
convenient way to circumvent this problem, and in addition it provides a new and unexpected 
perspective on the spectra of quantum graphs by allowing to compare their complexities. 

The expansion (5) is similar in spirit to the well-known EBK semiclassical quantization 
formula [4-6]. Given the quantum number n, Eq. (5) provides an individual expansion of 
the corresponding energy eigenvalue E^- In the same spirit EBK theory provides individual 
energy eigenvalues for a given set of quantum numbers by quantizing action integrals on 
tori. Thus the two methods are similar in the sense that both provide explicit values for 
the eigenenergies simply by plugging an integer (or a set of integers) into a known formula. 
This superficial similarity notwithstanding the underlying physics of the two methods is 
completely different. EBK relies on a simple, integrable structure of the underlying classical 
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dynamics based on (dynamical) symmetries whereas our method of exphcitly solving for the 
spectrum of quantum graphs relies on the construction of a network of spectral separators. 

The complexity of the expansion (5) compared to the EBK quantization formula reflects 
the complexity of the classical periodic orbit structure of quantum graphs. Moreover, the 
solution scheme shown above demonstrates that the spectral complexity of quantum graphs 
can be qualitatively different for different quantum graphs. According to this scheme, resolv- 
ing the irregular spectra may not amount to something as simple as redefining the expansion 
coefficients and the frequencies in Eq. (5). Hence, further generalization and simplification 
of the individual quantum eigenvalue quantization scheme outlined above will most likely 
prove to be highly nontrivial. Apparently, one encounters a whole hierarchy of complexities 
of the quantum spectra, even for such simple systems as the quasi one-dimensional quantum 
graphs. 

Concluding this section we would like to make a few comments on the comparison be- 
tween our analytical methods and standard numerical methods for computing eigenvalues 
of quantum graphs. In Ref. [11] we argued that there is an important conceptual difference 
between analytical and numerical methods. For instance analytical methods, such as ours, 
provide the solution of a whole class of objects simultaneously, whereas numerical methods 
address specific solutions of specific cases, one by one. In this sense analytical solutions are 
much more powerful than numerical solutions. In addition, our explicit analytical solutions 
of scaling quantum graphs are exact, whereas the accuracy of a numerical solution is bound 
by the word length of the computational device used, or, in case of "infinite-accuracy" al- 
gorithms, by the time one is willing to wait for the solution. Even if one is content with 
the numerical computation of finite accuracy, finite stretches of spectral eigenvalues, there 
are at least two situations, that require auxiliary analytical input: (i) the computation of 
spectral points for very large root number n and (ii) the computation of complete spectra. 
A discussion of both cases can be found in Ref. [11]. To this discussion we would hke to 
add the following recent development concerning the topic of complete spectra, which were 
required for the experimental and theoretical investigations of Ref. [47]. It was argued in 
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Ref. [47] that even a single missing state would have invalidated the experimental results 
reported in Ref. [47] . Certifying completeness of the experimental spectrum was only possi- 
ble with the help of numerical support, which itself used auxiliary analytical input to certify 
the completeness of the numerical spectra. This example illustrates clearly that the require- 
ment of complete spectra is not some idle academic pursuit, but that the need for complete 
spectra, and thus for analytical spectral methods, occurs in real-life situations, including 
experimental physics. 



In summary, we solved the spectral problem of scaling quantum graphs by deriving 
explicit, exact expressions for each individual energy eigenvalue E„ of the graph. On the 
level of the spectral equation our procedure for determining the energy eigenvalues also 
defines a method for solving analytically and explicitly a class of transcendental equations. 
This in itself is surprising and may have applications in pure mathematics, in particular in 
the theory of almost periodic functions [29] . 

The authors gratefully acknowledge financial support by NSF grant PHY-9984075. Work 
at UCSF was supported in part by the Sloan and Swartz Foundations. 

IX. APPENDIX A: PROOF OF "ONE ROOT PER ROOT CELL" 

Here we provide a proof for the statement (see Sec. II) that one and only one root A;„ of 
Eq. (19) is found in the root interval kn-i < k < kn, where kn are the root separators defined 
in Eq. (23). In order to simplify our task we scale and shift the argument k in Eq. (19), 



VIII. SUMMARY AND CONCLUSIONS 
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(/c + 7r7o). 
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and prove without loss of generality that 
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i=l 
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has precisely one zero kn in each interval /„ = {j^n-i, J^n), n & Z, i/n — "mr, if the regularity 
condition (22) is fulfilled. 
We start by showing that 



g(x) := 



1 - 



Dill «i COs{uJiX + Pi 



< 1 



for all X. The proof is straightforward. Defining 0^ = ujiX + Pi we have 



1 - 
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ai cos(0j 



> 1 



TV 



cos(6i 



> 1 - 



TV 



and, with Eq. (57), 



9{x) < 



Ej=i l«i sin(ei 
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(56) 



> 1 - > 0, (57) 



— r2<l- 
1 - [E»=i |aiCOs(ei)|J 

We now complete the proof in six steps. 

(i) We observe that \^{x)\ < E^Ii < 1 for all x. 

(n) We use (i) to show that the end points Un of /„ are not roots of Eq. (55): |-F(t', 

|(_l)n_$(^^)|>l_|$(^^)|>0. 

(iii) In In we define a new variable ^ according to 

X ^ Un + ^, < ^ < TT. 



(58) 



(59) 



Inserting Eq. (59) into Eq. (55) we see that in 7„ the spectral function F{x) is identical with 



/„(0 = (-!)" COs(0-¥'n(0, 



(60) 



where 
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N 

Vn{i) = XI cos(a;i^ + A + n'K'^i) . (61) 

i=l 

(iv) Because of (i) we have signF(i/„) — ("I)'*- We use this fact to show: 
signF(i/„)F(i/„+i) = (_i)2n+i _ _i Since F is continuous, this proves that there is at 
least one root of F in every /„. 

(v) According to (iii) and Eq. (60) the roots of F in satisfy (— l)'*cos(^) = (finiO^ 

^ = hn{0, (62) 

where /i„(^) = arccos[(— 1)"(^„(^)]. Therefore, roots of F are fixed points of 

(vi) In 7„, because of Eq. (56): 

- -| 2 

X)i=i Oio;! sin(a;i^ + A + rnrui) 

[h'niOr-^-^ ^ < 1- (63) 

1 - I]i=i «i cos(a;j^ + A + ri7ra;j) 

Prom Eq. (63) we obtain 

h'^iO < 1 in /„. (64) 

Because of Eq. (64) it now follows immediately that Eq. (62) has only a single fixed point. 
This is so since (iv) guarantees the existence of at least one fixed point ^* of Eq. (62). But 
because of Eq. (64) there cannot be any other, since Eq. (64) guarantees that |^ — /in (01 
increases monotonically to both sides of Consequently Eq. (62) has one and only one 
fixed point. Since, because of (v), the fixed points of Eq. (62) are the roots of Eq. (55) in 
/„, we showed that Eq. (55) has precisely one root in each root interval /„. 



X. APPENDIX B: CONVERGENCE OF PERIODIC ORBIT EXPANSIONS FOR 

INDIVIDUAL SPECTRAL POINTS 

Here we show that our explicit spectral formulas converge, and converge to the correct 
spectral eigenvalues. For the zeros of (19) we define the spectral staircase 
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N{k) = J2^ik-k), (65) 

i=l 

where 9{x) is Heavyside's 9 function (7). Based on the scattering quantization approach it 
was shown elsewhere [1] that 



1 °° 1 

N{k) = N{k) + -ImTr 7 (66) 

TT t 

1=1 

where 



— - (/^ + l + 7o), (67) 

TT 

and <S'(A;) is the unitary scattering matrix of the quantum graph. Since, according to our 
assumptions, S{k) is a finite, unitary matrix, existence and convergence of Eq. (66) is guar- 
anteed since in the eigenangle representation Eq. (66) involves nothing but the Fourier sums 
Ym^i sm{la{k))/l, which according to Ref. [26], formula 1.4411, converge to [tt — a{k)]/2 
mod 27r. Therefore, N{k) is well-defined for all k. Since S{k) can easily be constructed 
for any given quantum graph [1,10], Eq. (66) provides an explicit formula for the staircase 
function (65). This expression now enables us to explicitly compute the zeros of Eq. (19). 

In Appendix A we proved that exactly one zero kn of Eq. (19) is located in 7„ = {kn-i, kn)- 
Integrating A^(A;) from kn-i to kn and taking into account that N{k) jumps by one unit at 
k — kn, we obtain 

N{k) dk = N{kn-l)[kn - kn-l] + N {kn)[kn - kn]. (68) 

Solving for kn and using N{kn-i) — n — 1 and N{kn) — n, we obtain 

A;„ = ^ (2n + At + 7o) - [ N{k)dk. (69) 

"^0 Jkn-l 

Since we know N{k) explicitly, Eq. (69) allows us to compute every zero of Eq. (19) explicitly 
and individually for any choice of n. The representation (69) requires no further proof since, 
as mentioned above, N{k) is well-defined everywhere, and is integrable over any finite interval 
of A;. 
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Another useful representation of kn is obtained by substituting Eq. (66) with Eq. (67) 

into Eq. (69) and using — 7r[n + /x + 1/2 + 7o]/S'o: 

1 °° 1 

K^K ImTr / \"-S\k)dk. (70) 



1 oo 

kn = kn ImTrV- / S\k) dk. (71) 



In the eigenangle representation of the (S-matrix it is trivial to show by direct calculation 

that integration and summation can be interchanged in Eq. (70) and we arrive at 

I— I ' kn- 

In many cases the integral over S\k) can be performed explicitly, which yields explicit 
representations for 

Finally we discuss explicit representations of kn in terms of periodic orbits. Based on 
the product form of the S matrix [10] the trace of S\k) is of the form 

TrS{k)' = J2 Dn,nU,,,,,D^,,,U,,,, . . . D^^.^U,^,, = J] exp {zL^^\l]k} , (72) 

ji-jl meP[l] 

where P[l] is the index set of all possible periodic orbits of length / of the graph, Am[l] is 
the weight of orbit number m of length I, computable from the matrix elements of U, and 
Lto^[/] is the reduced action of periodic orbit number m of length I. Using this result we 
obtain the explicit periodic orbit formula for the spectrum in the form 



2 1 giim'Wfc^ 

kn = kn - -Im V- } Am[l] — z^r— sm 

1=1 m€P[l] ^'•^ 



(73) 



Since the derivation of Eq. (73) involves only a resummation of Tr S*' (which involves only a 
finite number of terms), the convergence properties of Eq. (71) are unaffected, and Eq. (73) 
converges. 

Reviewing our logic that took us from Eq. (69) to Eq. (73) it is important to stress that 
Eq. (73) converges to the correct result for kn- This is so because starting from Eq. (69), we 
arrive at Eq. (73) performing only allowed equivalence transformations. This is an important 
result. It means that even though Eq. (73) may only be conditionally convergent, it still 
converges to the correct result, provided the series is summed exactly as specified in Eq. (73). 

35 



The summation scheme specified in Eq. (73) means that periodic orbits have to be summed 
according to their symbohc lengths [39,40] and not, e.g., according to their action lengths. 
If this proviso is properly taken into account, Eq. (73) is an explicit, convergent periodic 
orbit representation for kn that converges to the exact value of kn- 
it is possible to re- write Eq. (73) into the more familiar form of summation over prime 
periodic orbits and their repetitions. Any periodic orbit m of length / in Eq. (73) consists 
of an irreducible, prime periodic orbit m-p of length l-p which is repeated i/ times, such that 

I = ul-p. (74) 

Of course v may be equal to 1 if orbit number m is already a prime periodic orbit. Let us 
now focus on the amplitude Am[l] in Eq. (71). If we denote by A^-p the amplitude of the 
prime periodic orbit, then 

A„,[l] = Ir A';,^. (75) 

This is so, because the prime periodic orbit mp is repeated u times, which by itself results in 
the amplitude A'^^ . The factor l-p is explained in the following way: because of the trace in 
Eq. (71), every vertex visited by the prime periodic orbit m-p contributes an amplitude A^^ 
to the total amplitude Am[^]- Since the prime periodic orbit is of length Ip, i.e. it visits l-p 
vertices, the total contribution is l-p A^^ . Finally, if we denote by Lm-p the reduced action 
of the prime periodic orbit m-p, then 

LtM--Ltl. (76) 
Collecting the results (74) - (76) and inserting them into Eq. (73) yields 



m-p J-imp i/=i 



rnp 



(77) 



where the summation is over all prime periodic orbits mp> of the graph and all their repeti- 
tions V. It is important to note here that the summation in Eq. (77) still has to be performed 
according to the symbolic lengths I — vl-p of the orbits. 
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In conclusion we note that our methods generahze and can be used to obtain any differ- 
entiable function directly and explicitly Integrating over f'{k)N{k) we obtain 

fikn) = nfCkn) -in- l)f{kn-i) - I fik) N{k) dk. (78) 

According to the same logic that led to Eq. (73), we obtain 

2 °° 1 

f{kn) = nfik) -{n- l)f{kn-i) - -Im J] y J] A^[l] Gn(i^^^[^]), (79) 

1=1 m.eP[l\ 

where 

Gn{x) = r f{k)e'^^dk. (80) 

This amounts to a resummation since one can also obtain the series for kn first, and then 
form f{kn). 
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